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Abstract 

We study the evolution of a flat Friedmann- Robertson- Walker Universe, filled with a bulk viscous 
g ; co.molog.ca! fluid, in the presence of tone varying "constants" . The dimensional analys.s of the 

model suggests a proportionality between the bulk viscous pressure of the dissipative fluid and the 
energy density. On using this assumption and with the choice of the standard equations of state 
\ for the bulk viscosity coefficient, temperature and relaxation time, the general solution of the field 

equations can be obtained, with all physical parameters having a power-law time dependence. The 
J> ■ symmetry analysis of this model, performed by using Lie group techniques, confirms the unicity of 

in . 

the solution for this functional form of the bulk viscous pressure. In order to find another possible 

O 

solution we relax the hypotheses assuming a concrete functional dependence for the "constants" . 
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I. INTRODUCTION 



Since the pioneering work of Dirac 1 , who proposed, motivated by the occurrence of large 
numbers in Universe, a theory with a time variable gravitational coupling constant G, cos- 
mological models with variable G and nonvanishing and variable cosmological term have 
been intensively investigated in the physical literature 2-14 . 

Dissipative thermodynamic processes in cosmology originating from a bulk viscosity are 
believed to play an important role in the dynamics and evolution of the Universe. Misner 15 
suggested that large-scale isotropy of the universe observed at the present epoch is due to 
the action of neutrino viscosity which was not negligible when the universe was less than a 
second old. There are a number of processes responsible for producing bulk viscosity in the 
early universe, such as the interaction between radiation and matter 16 , gravitational string 
production 17 " 18 , viscosity due to quark and gluon plasma, dark matter or particle creation 

24,33-34 

In order to study these phenomena, the theories of dissipation in Eckart-Landau 
formulation 19-20 , who made the first attempt at creating a relativistic theory of viscos- 
ity. However, these theories are now known to be pathological in several ways. Regardless 
of the choice of equation of state, all equilibrium states in these theories are unstable. In 
addition, as shown by Israel 21 , signals may be propagated through the fluid at velocities 
exceeding the speed of light in contradiction with the principle of causality. These problems 
arise due to the first order nature of the theory since it considers only first-order devia- 
tions from the equilibrium, leading to parabolic differential equations, and hence to infinite 
speeds of propagation for heat flow and viscosity. While such paradoxes appear particularly 
glaring in relativistic theory, infinite propagation speeds already constitutes a difficulty at 
the classical level, since one does not expect thermal disturbances to be carried faster than 
some (suitably defined) mean molecular speed. Conventional theory is thus applicable only 
to phenomena which are "quasi-stationary" i.e. slowly varying on space and time scales 
characterized by mean free path and mean collision time 21 . This is inadequate for many 
phenomena in high-energy astrophysics and relativistic cosmology involving steep gradients 
or rapid variations. These deficiencies can be traced to the fact that the conventional theories 
(both classical and relativistic) make overly restrictive hypothesis concerning the relation 
between the fluxes and densities of entropy, energy and particle number. 
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A relativistic second-order theory was found by Israel 21 and developed by Israel and 
Stewart 22 into what is called 'transient' or 'extended' irreversible thermodynamics. In this 
model, deviations from equilibrium (bulk stress, heat flow and shear stress) are treated as 
independent dynamical variables leading to a total of 14 dynamical fluid variables to be 
determined. The solutions of the full causal theory are well behaved for all times. Hence the 
advantages of the causal theories are the following 23 : 1) for stable fluid configurations, the 
dissipative signals propagate causally 2) unlike Eckart-type's theories, there is no generic 
short-wavelength secular instability in causal theories and 3) even for rotating fluids, the 
perturbations have a well-posed initial value problem. Therefore, the best currently available 
theory for analyzing dissipative processes in the Universe is the full Israel-Stewart causal 
thermodynamics. 

Due to the complicated nonlinear character of the evolution equations, very few exact 
cosmological solutions of the gravitational field equations are known in the framework of 
the full causal theory. For a homogeneous Universe filled with a full causal viscous fluid 
source obeying the relation £ ~ p% , exact general solutions of the field equations have been 
obtained earlier 24-33 . In this case the evolution of the bulk viscous cosmological model can 
be reduced to a Painleve-Ince type differential equation. It has also been proposed that 
causal bulk viscous thermodynamics can model on a phenomenological level matter creation 
in the early Universe 24-33 . 

In this paper, we consider the evolution of a causal bulk viscous fluid filled flat FRW 
type Universe, by assuming the standard equations of state for the bulk viscosity coefficient, 
temperature and relaxation time and in the presence of time varying constants (G, c and 
A). In order to obtain some very general properties of this cosmological model with variable 
constants, we shall adopt a method based on the studies of the symmetries of the field 
equations. As a first step, we shall study the field equations from dimensional point of 
view. The dimensional method provides general relations between physical quantities and 
allows us to make some definite assumptions on the behavior of thermodynamical quantities 
as well as on the equation of state for the bulk viscous parameter. In particular, we find 
that, under the assumption of the conservation of the total energy of the Universe, the bulk 
viscous pressure of the cosmological fluid must be proportional to the energy density of the 
matter component and for 7 = 1/2 (where 7 stands for the bulk viscous parameter), we 
find that G/c 2 must remain constant in spite of considering both "constants" as functions of 
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time t. On using this assumption, the gravitational field equations can be integrated exactly, 
leading to a general solution in which all thermodynamical quantities have a power-law time 
dependence. 

As we have been able to find a solution through dimensional analysis, it is possible that 
there are other symmetries of the model, since dimensional analysis is a reminiscent of scaling 
symmetries, which obviously are not the most general form of symmetries. Hence, we shall 
study the model through the method of Lie group symmetries, showing that under the 
assumed hypotheses of the proportionality of bulk viscous pressure to the energy density, 
there are no other solutions of the field equations. 

The paper is organized as follows: In section (II), we outline the equations of the model, 
define all the quantities as well as their equations of state and fix the notation. In section 
(III), we study our equations from the dimensional point of view (subsection III A). This 
study allows us to make two simplifying hypotheses which we use to integrate the equations. 
In subsection (III B) , we obtain the equation of state for the bulk viscosity through the 
dynamical approach. Section (IV) presents a naive method to study our model taking into 
account the previous hypotheses. This technique allows us to obtain a complete set of 
solutions as well as to arrive at some interesting conclusions. In section (V), we check the 
solutions obtained in the previous sections are correct and show that direct integration does 
not give more information about the solutions. As Dimensional Analysis is a manifestation 
of symmetries, in section (VI) we show that under the imposed hypothesis there are no 
more solutions for the field equations than ones already obtained. In section (VII), we try 
to relax the hypotheses in order to seek other solutions but (un)fortunately we only recover 
our previous solutions. We end the paper by summarising our findings in section (VIII). 

II. THE MODEL 

Following Maartens 35 , we consider a Friedmann- Robertson- Walker (FRW) Universe with 
a line element 



filled with a bulk viscous cosmological fluid with the following energy-momentum tensor: 



ds 2 = c 2 dt 2 - f 2 {t) (dx 2 + dy 2 + dz 2 ) , 



(1) 



T? = (j> + p + IL) Ui u k -{ V + U)5l 



(2) 
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where p is the energy density, p the thermodynamic pressure, II is the bulk viscous pressure 
and Ui is the four velocity satisfying the condition u(u % = 1. The number 4-flux and the 
entropy 4-flux take the form: 

N* = nu\ (3) 
& = sN { - ( j u\ (4) 



where n is the number density, s the specific entropy, T > the temperature, £ the bulk 
viscosity coefficient and r > the relaxation coefficient for transient bulk viscous effect 
(i.e. the relaxation time). The fundamental thermodynamic tensors (2-4) are subject to 
the dynamical laws of energy-momentum conservation, number conservation and the Gibb's 
equation: 

T* k = 0, (5) 

n; = o, (6) 

Tds = d(^) +pd(-) . (7) 



The equations (2-4) and (5-7) imply: 

rsj = -n f 3^ + + ^rn , (8) 

by equations (5-7) and (8), the simplest way (linear in II) to satisfy the if -theorem (i.e. 
for the entropy production to be non-negative, = ^ > 22 > 36 ~ 37 ) leads to the causal 
evolution equation for bulk viscosity given by 35 

(9) 

In eq.(9), e = gives the truncated theory (the truncated theory implies a drastic 
condition on the temperature), while e = 1 gives the full theory. The non-causal theory has 
r = 0. 

The growth of the total commoving entropy £ over a proper time interval (to,t) is given 
by 35 : 

E(f) - E (to) = ~ t^^dt, (10) 



k B J t T 
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where kg is the Boltzmann's constant. The Einstein gravitational field equations with 
variable G, c and A are: 

1 „ 87rGm 

Rik - ^9ikR = -^rnr T * + A (*)^<*- ( n ) 

Applying the covariance divergence to the second member of equation (11) we get: 

div [^Ti + 5{AJ = 0, (12) 

/4c, G A _,• c 4 ^ A 7 - , . 

that simplifies to: 



c G / * 8ttG 



Ac 4 G c 

p + 3 (p + p) H + 3HU = -— - p- - 4p-, (14) 

OTVLr Cr C 

where if stands for the Hubble parameter (H = ///). The last equation may be written in 
the following form: 

Ac 4 C r 

p + 3(p + p)H + 3HU+—- + p--4p-=0. (15) 

07TG Lr c 

Therefore, our model (with FRW symmetries) is described by the following equations: 



8ttG 



2H + 3H 2 = — — (p + IT) + Ac 2 , (16) 



3H 2 = -— p + Ac 2 , (17) 
c z 

Ac 4 G c 

p + 3{p + p + Il )H = - - p +4p-, (18) 

07T(jr Ct c 

rli + n = - |rn f 3fT + I - | - (19) 

In order to close the system of equations (16-19) we have to give the equation of state 
for p and specify T, £ and r. As usual, we assume the following phenomenological (ad hoc) 
laws 35 : 

P = up, (20) 

t = fc 7 P 7 , (21) 

T = D 5 p 5 , (22) 

T = frr 1 = ky~\ (23) 



where < to < 1, and A; 7 > 0, -D,? > are dimensional constants, 7 > and 5 > are 
numerical constants. Eqs. (20) are standard in cosmological models whereas the equation 
for t is a simple procedure to ensure that the speed of viscous pulses does not exceed the 
speed of light. These are without sufficient thermodynamical motivation, but, in absence 
of better alternatives, we use these equations and expect that they will at least provide an 
indication of the range of possibilities. For the temperature law, we take, T = D$p 5 , which 
is the simplest law guaranteeing positive heat capacity. 

In the context of irreversible thermodynamics p, p, T and the particle number density n 
are equilibrium magnitudes which are related by equations of state of the form p = p(T, n) 
and p = p(T,n). From the requirement that the entropy is a state function, we obtain the 
equation 



For the equations of state (20-23) this relation imposes the constraint 5 = so that 
< 5 < 1/2 for < oj < 1, a range of values which is usually considered in the physical 
literature 38 . 

The Israel-Stewart-Hiscock theory is derived under the assumption that the thermody- 
namical state of the fluid is close to equilibrium, that is, the non-equilibrium bulk viscous 
pressure should be small when compared to the local equilibrium pressure |n| << p = up 
39-40 _ jf this condition is violated then one is effectively assuming that the linear theory 
holds also in the nonlinear regime far from equilibrium. For a fluid description of the matter, 
this condition should be satisfied. 

Therefore, with all these assumptions and taking into account the conservation principle, 
i.e., div(T-) = 0, the resulting field equations are as follows: 




(24) 



2H + 3H 2 



3H 2 




(26) 
(27) 



(25) 




Ac 4 



+ ^-4p- = 0, 



(28) 



8ttG 




(29) 
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where W = (^Jf). 

III. HYPOTHESES AND EQUATIONS OF STATE 

In this section, we study the field equations in a dimensional way as well as with a dynam- 
ical systems approach in order to obtain relationships between various physical quantities, 
simplify the field equations and to obtain an adequate equation of state for the bulk viscous 
parameter. 

We use dimensional considerations to re-write the field equations in a dimensionless way 
and show that some interesting relations can be obtained as a result. We find that the 
bulk viscous pressure has the same behaviour as the energy density and if 7 7^ 1/2, the 
relationship G/c 2 should vary if we demand that our equations remain gauge and scale 
invariant. We find that with the special case 7=1/2 the relationship G/c 2 = const, but 
in such a way that both "constants" vary. Dimensional analysis suggests that 7 = 1/2 is a 
very special value for the bulk viscous parameter. Using the dynamical systems approach, 
we shall show that for this value of 7, our model approximates the dynamics of a perfect 
fluid and for large time, the dynamics of the Universe follows that of a flat FRW model. 

A. Dimensional considerations 

The 7r — monomia is the main object in dimensional analysis. It may be defined as the 
product of quantities which are invariant under change of fundamental units, ir — monomia 
are dimensionless quantities, their dimensions are equal to unity. The dimensional analysis 
has structure of Lie group 41-44 . The it — monomia are invariants under the action of the 
similarity group. We must mention here that the similarity group is only a special class of 
the group of all symmetries that can be obtained using the Lie method. For this reason, 
when one uses dimensional analysis only one of the several possible solutions to the problem 
is obtained. 

The equations (25-29) and the equation of state (20-23) can be expressed in a dimension- 
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monomia: 








Gpt 2 




Gilt 2 




Gpt 2 


c 2 


c 2 


c 2 


TT 
11 

V 


7T 5 


. J_ 
lit 


7T6 = 


- T - = tH 
t 


t 




? 


7T 9 = 


T 
1 




rp 


D s p s 


P 
P 




= kc 2 t 2 




Ac 

Gp 



n ± = ^ = — — n 3 = —g- (30) 

(31) 

717 =k^ ** = T P n9= D^ (32) 

7Ti = - 7Tn = AcV 7Ti2 = 7^- (33) 

The following relations can be obtained from the n — monomia. 

1. From 71"!, 7r 2 , 7r 3 , 7r 4 and 7r 10 , we see that 

p oc p oc II, (34) 

and thus, [II] = [p] = [p]. Therefore, these physical quantities have the same dimen- 
sional equation and exhibit the same behaviour in order of magnitude. 

2. On using 7r 7 and 7r 8 , we obtain the following relation 

n 8 = ^^^ P =(k-H) b (35) 

and from 7r 3 and 7r 8 

Gpt 2 _ k lP ^ c 2 ft V /(7 ^ 1} _ ^ _ fW (36) 



t Gt 2 \ /c^ y G~ 



where b = 1/(7 — 1). When 7 = 1/2, we obtain the relationship k 2 = c /G, which 
indicates that these relations remain constant, i.e., G and c can vary but in such a 
way that c 2 /G remains constant. If 7 ^ 1/2, the "constants" G and c must vary if we 
want that our equations to remain gauge invariant or scale invariant. In a reasonable 
physical approach, we need to impose the condition G/c 2 = constant, and in the 
process, uniquely fixing the value of the coefficient 7. The dimensional analysis not 
only fixes the value of 7, but also predicts the value of the proportionality coefficient 
as being k 2 = c 2 /G (in natural units k 2 = 8tt). 
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Therefore, the solutions that Dimensional Analysis suggests us are the following: 



P oc (k~H) , (37) 

p = up on (k~ l t) b , (38) 

n oc (k~H) b = np, (39) 

A oc c" 2 t- 2 , (40) 

% oc k^ h -\ (41) 
where k is a numerical constant. From the equation of state, we obtain 

£ = k^p* oc A; 7 (A^t) 67 , (42) 

T = £> 5 / oc fc 1 *) 56 with 5 = (43) 
v ' ' uj + 1 

r = ^p- 1 = ^ 7 (A^ 1 *)^ , i. e. r = t (44) 



So far, we have seen that with this simple method we can obtain the behaviour of all the 
physical quantities. We would like to emphasize that we only have compared tt — monomia 
and used the equation of state to obtain these results. In section IV, we will present a naive 
way to work with another dimensional method. 

B. Structural stability of bulk viscous cosmological models 

In this subsection, we consider the analysis of bulk viscous cosmological models using 
dynamical systems theory. Dynamical systems have already been used in the study of 
causal viscous fluids 28-30 , non-causal viscous fluids 45 , G2 cosmologies 46-47 and magnetic 
fields in scalar field cosmology 48 (for a review of dynamical systems in cosmology, see 49_50 ). 
In this method, the governing equations of the model are considered as a finite system of 
autonomous ordinary differential equations. By using this method, we can determine the 
arbitrary exponent in the state equation of the bulk viscosity coefficient. 

As a basic physical requirement, we assume that the causal bulk viscous model approx- 
imates the dynamics of a perfect fluid and for large time, the dynamics of the Universe 
follows that of a flat FRW model. Our assumption is supported by the fact that the viscous 
pressure will decay faster than the thermodynamic pressure 39 and hence the behaviour of 
our model is equivalent to that of a model with a perfect fluid. On a large scale, the Uni- 
verse is well described by the classical flat FRW models. It is natural to assume here that 
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perturbations introduced by a dissipative parameter like bulk viscosity does not modify the 
cosmological evolution even at early times. Viscous dissipative terms can significantly mod- 
ify the thermodynamical behaviour of the early cosmological fluid, but, we do not expect a 
major modification in the scale factor. With these assumptions, we can use mathematical 
techniques of dynamical systems theory to select a correct equation of state for the bulk 
viscosity coefficient, i.e. the values of the parameter 7, in such a way that the viscous model 
describes a homeomorphic dynamics of the flat FRW model. A natural environment to study 
the first order autonomous ordinary differential equations 

^ = P(x, y,z),^ = Q(x, y,z),-^ = R(x, y, z), (45) 

is a three-dimensional differential manifold M to every point p = (x, y, z) of which the 
tangent space T P M is associated. The vector field X G X(M) / X(p) = (P(p), Q(p), R(p)) G 
T p M of class C 1 . The set X(M) of all such vector fields on M is called the space of differential 
equations on M. The solution curves, <pt{p) (also called phase trajectories) of a vector field 
X(p) on M define a one-parameter group of transformations <p t '■ M — * M where t G (a, b) . 
If M is compact, then i e M and 4> t is called the dynamical system. 

Two vector fields X, Y e X(M) are said to be topologically equivalent if there is a 
homeomorphism h : M — > M preserving the orientation of the solution curves of X into 
those of Y. Let D C M be a compact subset of M. We say that X G X(D) is structurally 
stable if there is a neighborhood N(X) of X in X(D) such that for every Y C N(X), Y 
is topologically equivalent to X, i.e. there is an orientation preserving homeomorphism 
transforming the phase trajectories of X into those of Y. 

In the following, we seek spatially homogeneous and isotropic bulk viscous solutions of the 
gravitational field equations which, in a structurally stable way, approximate the dynamics 
of ordinary (perfect fluid) FRW models. More precisely, we will look for such spatially 
and isotropic solutions which, after being perturbed by the dissipative parameter, yield the 
dynamics topologically equivalent to those of FRW models. If we take into account the state 
equations (20-23), the gravitational field and the bulk viscous pressure evolution equations 
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can be rewritten in the following form: 



H = -H 



,_(3(« + l)-2) 1 1 
6 ^23 



P 

ri 



-3i/((cu + i)p + n) 



" r + 3p# + ^n 



3#- W 



with 



Solving 



H 2 = -p h -A. 

3 F 6 3 



if = 0, p = 0, n = o, 
for H, p and n, we obtain from (47) 

n = -(w + i)p, 

and taking into account this result, on using (46) we arrive at: 



H 2 = -(p + A) 



(46) 
(47) 
(48) 

(49) 
(50) 



(51) 



With the use of Eq. (51), the energy density p can be obtained from the equation (note 
that p = 0) 



n 



' r + 3 P # + in 



3H- W P - 

P. 



0. 



and is given by 



or alternately 



p{k~ 2 (u; + l) 2 P 1 -^-(3)=(3A J 



where j3 = 3 ^1 — ^ ^ ^ . As we can see 



, ,2 P 

9 ~ 1 ^ 



7 



(52) 



(53) 



(54) 



(55) 



since < h Vo;. For this value, we have p = A = 0. 

Therefore, as has already been pointed out earlier 24-26 , causal bulk viscous cosmo logical 
models with 7 = 1/2 describe a two-phase evolution of the Universe. The Universe is 
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born from a zero energy density vacuum state and in the first phase the energy density is 
increasing to a maximum value. This period corresponds to a phase of matter creation. 
After reaching a maximum value, due to the expansion of the Universe, the energy density 
becomes a monotonically decreasing function of time and the second phase describes a 
standard evolution of the bulk viscous causal cosmological fluid. Hence, in this framework, 
bulk viscous processes can model matter creation in the early Universe. The ultra-stiff case 
of the Zel'dovich fluid, to — 1 also gives p = 0, but from this value we cannot find the value 
for 7. 

The structurally stable approximation of the system perturbed by the viscous parameter 
must have the same critical points as the unperturbed system. The unperturbed system has 
no critical points except the point (H,p) = (0,0). Therefore it is observed that the only 
structurally stable approximations to the flat FRW solutions are those with £ = k^p 1 ^ 2 , i.e. 
7 = 1/2. 

As already have been pointed out by M. Szydlowski and M. Heller 51 in the study of 
non-causal viscous fluids it is possible to use a Lie group technique in order to determine 
an adequate equation of state, in particular theses authors find that there is a only one 
symmetry for the field equations if 7 = 1/2. 

IV. A NAIVE METHOD 

In this section, we use dimensional analysis to obtain a complete set of solutions for the 
field equations with the assumptions div{T) = and 11 = xp and x e Mr. This leads to 

P = A^r 3 ^ +1+ ^ or p = A^f~ a , (56) 

valid for V7. The set of governing parameters are SDT = ^{A^ >3< , k 7 , t} , where [A^] = 
L 3(w+i+^)-i MT -2 ? = L -i/2 M i/2 and ^ = T Uging thege5 we bt a in the following 
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relations: 



G oc A^k^r 4 -^), (57) 

c oc AS^k^T 1 -^, (58) 

p oc fc~ 6 ~V _1 oc n, (59) 

£ = fc 7 p 7 > r = ep-" 1 (60) 



/ oc AS£fcT r 1 ™, (61) 



k B 9 oc A#Ffc? Q+1) ^-5^+1) + ^, (62) 

a^s oc ^Jfel^^^r^^r), (63) 

A oc AS^k^t^, (64) 

g = -ba - b - 1 (65) 

where, a = 3(u> + 1 + x) — 1, 6 = 7 — 1, sis the entropy and g is the deceleration parameter. 
Note that the current method is very different from the method employed in the previous 
section. 

It is observed that 

— = const. <^ b = ~ 7 = i (66) 
and for 7 = ~, we obtain the following results: 

G oc t" 2 (Sr), coct"(Sr), p oc i -2 oc IT, (67) 

/ oc t^ 1 ), oc t" 2 (Si), a" 1/4 s oc jK&Sj), (68) 

-4 a — 1 . 

A oc g = — — . (69) 

Thus, we see that (for 7 = |) G, c and A are decreasing functions on time despite the fact 
that the relation G/c 2 remains constant. In a previous work 52 , where we only consider G and 
A as time functions, we found that for 7 = 1/2, G behaves as a true constant i.e. G = const. 
and that A vanishes. In this new scenario, we find that G and c varies with time but in 
such a way that the relation G/c 2 = const, and A is a decreasing funtion on time. This 
consequence has some interesting implications in light of recent supernovae observations 
which suggest a small but non-zero value for the cosmo logical constant 53 . 

We would like to point out that these results are very similar to those in 54 where a per- 
fect fluid cosmological model with time- varying constants in the presence of adiabatic matter 
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creation was considered. This striking similarity of our results is not surprising as it has 
been suggested by Zeldovich 55 and later by Murphy 56 and Hu 57 that the introduction of vis- 
cosity in the cosmological fluid is a phenomenological description of the effect of creation of 
particles by the non-stationary gravitational field of the expanding cosmos. A non-vanishing 
particle production rate is equivalent to a bulk viscous pressure in the cosmological fluid or 
from a quantum point of view, with a viscosity of the vacuum. This is due to the simple cir- 
cumstance that any source term in the energy balance of a relativistic fluid may be formally 
rewritten in terms of an effective bulk viscosity. Zimdahl and others 40 have considered in 
detail the possibility that the bulk viscous pressure of the full Israel-Stewart-Hiscock theory 
may also be interpreted as an effective description for particle production processes. The 
particle creation process leads to considerable changes in the thermodynamical behavior of 
the Universe. If the chemical potential of the newly created particles is zero, fi = 0, then the 
non-vanishing bulk pressure II associated with an increase in the number of fluid particles 
satisfies formally the same equation as in the case of presence of a real dissipative bulk 
viscosity 24-33 . 

This coincidence can be easily explained by observing the set of the governing parameters. 

In the case of a viscous fluid, this set is characterized by DJl = 9Jt {A^^, k 7 , t} , where [A UJ> J = 
i 3Ki+x)-i MT -2 ) = L -i/2 M i/2 and = T? whUe in the cage of adiabatic ma tter 

creation, this set is characterized by Wl = m{A ult p, B,t} with [A^^] = £3(«+i)(i-/J)-i MT~ 2 , 
[G/c 2 = B] = LM~ X and [t] = T, see 41 for the details. We can observe the extreme similarity 
between A^,* and Au$ and the coincidence between B and k~ 2 , which has been shown in 
the previous section as k 2 = c 2 /G. 

Another similarity of the present results with those obtained in 54 is that for adiabatic 
matter creation, we cannot impose such a mechanism with uj = on the equation of state 
for matter predominance. In the same way, we observe from eq. (62) that if we fix uo = 
into the equation of state the fluid temperature increases indefinitely with time. Thus, we 
arrive at a similar conclusion that, in the matter predominance era, we cannot consider the 
effects of the bulk viscosity in our fluid. This fact is in agreement with a previous result 
obtained in 58 where we studied a full causal bulk viscous model through the renormalization 
group method arriving to the conclusion that in the long-time asymptotics the cosmological 
model tends to a flat ideal (non-viscous) Friedmann type geometry. 

We also notice that the x— parameter, i.e. the causal bulk viscous effect, weakly perturbs 
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the FRW perfect fluid solution in such a way that the comoving entropy varies with t, 

while in the perfect fluid case, i.e. x = 0, the comoving entropy is constant. The same 

i 

considerations could be taken into account for the scale factor f(t) = f t 3< - 1 -^^+ 1+>! \ which 
is perturbed weakly for the bulk viscous parameter as we can see in the special case of 
7 = 1/2. The solutions obtained for the scale factor, energy density, entropy, etc. are similar 
"but different" to those obtained for a perfect fluid model where the x— parameter vanishes. 
We would like to point out that our model is thermodynamically consistent for the usual 
matter equations of state and valid for all 7— parameter, i.e. V7 G (0, 1) . In the same way, 
we can see that the viscous parameter helps us to get rid of the so-called entropy problem 
since in this model entropy varies with t. 

The bulk viscous pressure and the energy density of the cosmological fluid are proportional 
to each other and hence, the general evolution of II is qualitatively similar to the evolution 
of the thermodynamic pressure p, both obeying a similar equation of state. In fact, by 
defining an effective coefficient u e ff = uo — x, the equation of state of the cosmological 
fluid can be formally written as p e ff = p + IT = {w e ff — 1) p. However, since the bulk 
viscous pressure of the cosmological fluid must also satisfy the evolution equation (29), the 
resulting time evolution depends not only on w e //, but also, via the coefficients 5 and 7, on 
the equations of state of the bulk viscosity coefficient, £ = £ (p) and of the temperature, 
T = T (p). Therefore, even that formally one can introduce an effective pressure (including 
both p and II), obeying a cj-law equation of state, due to the extra-constraints imposed by 
the requirements of the causal bulk evolution, the general dynamics of the present model 
cannot be reduced to a perfect fluid model evolution. 

V. DIRECT INTEGRATION OF THE FIELD EQUATIONS 

As we have seen in the previous section, dimensional analysis suggests us that II oc p. We 
formalize this relationship in the following way, II = xp, and under physical considerations 
we take x G M~ and |x| <C 1 40 - In this section, we integrate the field equations with 
these assumptions and taking into account the conservation principle. This leads us to the 
following relationship 

p + 3 (u + 1 + x) pH = 0. (70) 
This trivially lead us to the well known relationship between the energy density p and 
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the scale factor / 

p = 4,r 3(w+1+x) or p = AJ~ a , (71) 

where a = 3 (a; + 1 + x). Now, taking into account the eq. (19) and simplifying it, we 
obtain, 

*P + ^ = -p-±*p(---w), (72) 
fc 7 p 7 1 a 2 \ a J 

thereby obtaining p = p(t). On simplifying further, we obtain, 

4; = T =► P = ^ 

where 



p = dk~H\ (73) 



^ = d = ( 7 _ l) K b and b = -. (74) 

From equation (71) we obtain: 

/=[-^fc?<~ 6 J , i-e. / oct^+i+xK7-D. (75) 

An important observational quantity is the deceleration parameter q = 4i{jj) — 1- The 
sign of the deceleration parameter indicates whether the model inflates or not. The positive 
sign of q corresponds to "standard" decelerating models whereas the negative sign indicates 
inflation. In our model, the deceleration parameter behaves as: 

?=-!-£ (76) 

Thus, for the parameters 7, u and x, the deceleration parameter indicates an inflationary 
behaviour when a/b > and a non- inflationary behaviour for a/b < —1. We proceed with 
the calculation of the other physical quantities as under: 

£ = oc A; 7 (dk- b t b Y = d^k~ b ^ b , (77) 
T = D 5 p s = D 5 (dk~ b t b ) S with 5 = (78) 
r = ^p" 1 = (dk~ b t b ) (l ~ l) , i.e. r = d^H (79) 

We see from r = d 7 " 1 ^ that this result is in agreement with the theoretical result obtained 
in 35 . For viscous expansion to be non-thermalizing, we should have r < t, or otherwise the 
basic interaction rate for viscous effects should be sufficiently rapid to restore the equilibrium 
as the fluid expands. The comoving entropy is 



m - s M = -i^i"^<-^«> IPO-*-"]- • (so) 

d~+ 5 k B D & L J to 
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We notice that the parameter x weakly perturbs the perfect fluid FRW Universe. When 
x = 0, the comoving entropy assumes a constant value and we recover the perfect fluid case. 
Finally, we will use the equations 



3H 2 = ^fp + Ac 2 , (81) 



Ac 4 G c 

--4- = 0, (82) 



8nGp G c 

to obtain the behaviour of the "constants" G, c and A. From (81), we obtain A and using it 
in eq. (82) we obtain 

3cc „, c 2 3HH p , . 

H + — = (83) 



AnxGp AnGp p 

Considering the previous results 



H = (3r 1 and p = dk~ h t\ (84) 



where (3 = —b/a. 

These results are used in eq. (83) to obtain the next ODE 



2 „2 + -2-b 



3/? 2 cH 



And Gk~ b 



c 1 



+ bj = 0. (85) 



We can see that the above equation verifies the relationship (a particular solution) 

k b c 2 

gt b+2 therefore G = gk b c 2 t- b - 2 , (86) 



g being a numerical constant. It is observed that if b = —2 i.e. 7 = 1/2, we obtain G/c 2 = B. 
We use this relation in eq.(85) obtaining 

where K c = and therefore, 

c = K c t K , (88) 

where K c is an integration constant and k = ^—1 — -^-^ . 

If we impose a new assumption (which is supported by observational considerations of A 
decaying with time and assuming a small but non-zero value) as 
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we can see from eq. (81) that 

'3(3 2 - I 



G 

72 



k b t -b-2 = Bk b t -b-2 (g 0) 



8ird 

Taking all these relations into consideration, we use in eq.(82) to obtain 

'(° + ±)-2°- b -±*=0, (91) 

AttB \c t J c t y J 

whose trivial solution is 

c = R c r», (92) 

where $ c is an integration constant and p = (j^^^^j ■ 

We have solved our model under the assumptions: div{T) = 0, n = xp, with x e M~ 
and valid for V7. It is observed that if 7 = 1/2 we obtain G/c 2 = const, as obtained earlier. 
Thus with the imposed hypotheses, we have obtained similar results. 

VI. LIE METHOD 

In this section, we study the field equations using the symmetry method. As we have 
discussed earlier, dimensional analysis is just a manifestation of scaling symmetry. However, 
this type of symmetry is not the most general form of symmetries 59 " 61 . Therefore, by 
studying the form of G(t) and c(t) for which the equations admit symmetries, we expect 
to uncover new integrable models. With these assumptions, we shall see that the solutions 
obtained in previous sections are recovered. 

We start with the assumption II = xp, with x £ M~ The bulk viscosity evolution equation 
can then be rewritten in the alternative form 

5^ + fc-V 1 " 7 = 30H, (93) 
P 7 

where /3 = (J - \) and 8 = (l - f ) . 

Taking the derivative with respect to the time of this equation and with the use of the 
next equation 

H = -4™§|p, (94) 
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obtained from the field equations (where a = (1 + u + x)), we obtain the following second 
order differential equation describing the time variation of the density of the cosmological 
fluid: 

P=P--Ap s p + B < ^§-p\ (95) 
P c 2 (t) 

where A = k ^^\ s = (1 - 7 ) an d B = ^f£. 

Equation (95) is of the general form. 

P = ^(t,p,p), (96) 

where i/>(t, p,p) = £- Ap s p + B^p 2 . 

We are going now to apply all the standard procedure of Lie group analysis to this 
equation (see 59 for details and notation) 

A vector field X 

X = at,p)d t + r){t,p)d p , (97) 

is a symmetry of (95) if 

- - # P + rju + (2%, - £«) p + (v P p - 26 P ) p 2 ~ £ PP P 3 + 



+ ( Vp - 2& - 3p£p) ^ - + (^ - 6) P - P 2 Q ip P = 0. (98) 

By expanding and separating (98) with respect to powers of p we obtain the overdeter- 
mined system: 

Z P p + p-% = 0, (99) 
V p- 2 - ^p" 1 + Vpp - 2C tp + 2Ai p p l ~ s = 0, (100) 
2rkp - Ztt + MtP s ~ 3B^p 2 - 27]^ + Asrip*- 1 = 0, (101) 



■ B £ I % ~ 2G ^ ) P 2 - 2B ^P + + (V P - 26) B^p 2 + A^p 1 "* = 0, (102) 



Solving (99-102), we find that 

t(p,t) = -ast + b, rj(p,t) = ap, (103) 

with the constraint 

G c (l-2s)a , 
G c o — ast 
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a and b being numerical constants. 

Thus, we have found all the possible forms of G and c for which eq. (95) admits sym- 
metries. There are two cases with respect to the values of the constant a; a = which 
correspond to G = const, c = const, and a^O which correspond to 

^- = K[ast + b]- 2+ K (105) 
cr 

with K a constant of integration. If s = |, which corresponds to 7 — |, then G and c remain 
constants. For this form of G and c eq. (95) admits a single symmetry 

X = (ast - b) d t - (ap) d p . (106) 

The knowledge of one symmetry X might suggest the form of a particular solution as an 
invariant of the operator X, i.e. the solution of 

#r = -r^v ( 107 ) 

This particular solution is known as an invariant solution (generalization of similarity 
solution). In this case 

p = Pot~^ (108) 

where for simplicity we have taken 6 = 0, and po is a constant of integration (note that 
s = (1 — 7) being 7 the bulk viscous parameter). 

We can apply a pedestrian method to try to obtain the same results. In this way, 
taking into account dimensional considerations, from the eq. (95) we obtain the following 
relationships between density, time and gravitational constant: 



G 

kZ X 8p't ^ 1, B—pt 2 ~ 1. (109) 

7 cr 



This last relationship is also known as the relation for inertia obtained by Sciama. From 
these relationships, we obtain 

p w B~ X K- X c 2 [ast + b]^ t~ 2 w B~ x K~ x c 2 [ast]~- . (110) 

We can see that this result verifies the relation 

fcrV*-i- (in) 
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Once we have obtained p, we can obtain / (the scale factor) from 

P = A^r 3(M ^ / = (A^H)^^ , (112) 
In this way, we find H and from eq. (26) the behaviour of A is obtained as: 



(3k 2 - 87rif(as)- 2+ ^p ) \ 



A = ^ S5 < 113 » 

where, k = 3^+1+. \}s • ^ n re pl ac i n g & U these results into eq. (28), we obtain the exact 
behaviour for c, i.e. 

I 1 - 2s\ 1 / I \ c , , 

+ 7 = — + 2 - 114 

c = c t- £ , (115) 



therefore 



where, e = y , L and A = — ~— . 

Therefore, the assumption II = up together with the equation (29) is very restrictive. 
Hence, we conclude that under these assumptions, the field equations do not admit any 
other solution. 

VII. RELAXING THE HYPOTHESES 

In the previous section, we have seen that our assumption II = xp together with the 
equation (29) does not give rise to any new solutions. This motivates us to investigate our 
assumptions further and relax the hypotheses we have used in obtaining solutions to the 
field equations. With the following assumptions obtained (suggested) from the sub-section 
III A, we would like to obtain an equation that helps us to determine the behaviour of the 
"constants" . The new hypotheses are: 



1. Relationship between G and c: 



2. A follows the law: 



% = Kk b 1 H 2+b . (116) 



/ FF 2 
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where, d 



3-1 



where K, I £ R are numerical constants and b = 1/(7 — 1). 

In this way, we allow the scale factor / to be an arbitrary function of time and it does 
not follow a power law i.e. / oc t a as in previous sections. From the eq. (81) and taking 
into account the assumptions, we obtain 

p = dH~\ (118) 
. Therefore from (27), we obtain 

n= — H' 2 ~ b H -d(uj + l)dH- b (119) 
3 

If we introduce these results into the equation (29), then 

H + KiHH + K 2 H~ l H 2 + K 3 H 3 = (120) 

where W = (^) and 

K 1= (3 + K ^), tf 2 = (-6-2 + W6) ,a-, = 3 (£±^- ^> + 1 ) ), 

(121) 

In order to solve eq. (120), we use the Lie method which results in the following overdeter- 
mined system: 

-Zhh + ZhKzH- 1 =0 (122) 
- V K 2 H- 2 + t]hh - 2£w + VHK2H- 1 + 2£ H K X H = (123) 
rjK, + 2 VtH - Ctt + ZtK x H + 3£ H K 3 H 3 + 2r /t A' 2 #- 1 = (124) 
3r]K 3 H 2 + rj tt - r] H K 3 H 3 + 2&K 3 H 3 + rj^H = (125) 

Solving (122-125), we find that 

f (t, if) = -at + 6, r](t, H) = aH (126) 

a and being numerical constants. For this form of £ and 77, eq. (120) admits a single 
symmetry 

X = (at - b)d t - (aH) d H . (127) 
Following the standard procedure, we find 

dt dH H = ^-, (128) 



£ (t, H) 7] (t, #) at - 6 
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which imply that 

f = (at-b)' (129) 

Thus, we obtain a power law solution as in previous section. To complete the set of equations, 
we have to take into account the next equation 

and replace it into all our results thereby leading to the next relationship, 

c = KH a (131) 

where a = l+ ^2A^ an< ^ ^ = ^ n Kd. In this way, we can see that 

G = KH 2(a+1)+b and A = lH 2{l ~ a \ (132) 

On using a straightforward calculation, we see that we have obtained the same solution 
even with different assumptions. The hypotheses considered earlier were relaxed in the 
above analysis with a view of obtaining a new set of solutions. This suggests that solutions 
obtained earlier are of a fairly general nature and do not really depend on any specific 
assumptions. 



VIII. CONCLUSIONS. 



In the present paper, we have studied a causal bulk viscous cosmological model with time 
varying constants. The dimensional analysis of the model allows us to make the following 
assumptions: the bulk viscous pressure is proportional to the energy density of the cosmo- 
logical fluid and that for 7 = 1/2, the relation G/c 2 remains constant in spite of considering 
"constants" G and c as functions of t. In addition, we have shown that for 7 = 1/2, the 
dynamics of our model is similar to that a perfect fluid FRW model. 

With the help of these assumptions, the general solution of the gravitational equations can 
be obtained in an exact form, leading to a power-law behavour of the physical parameters 
on the cosmological time. These solutions show us that the "constants" G, c and A are 
decreasing functions of time and that for 7 = 1/2, we obtain a model similar to a perfect fluid 
cosmological model with time- varying constants in the presence of adiabatic matter creation. 
This striking similarity of our results is not surprising as it has been pointed out by many 
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authors that the introduction of viscosity in the cosmological fluid is a phenomenological 

description of the effect of creation of particles by the non- stationary gravitational field of the 

expanding cosmos. We have shown that we cannot consider the bulk viscosity in the matter 

predominance era since the fluid temperature increases indefinitely with time. We arrive at 

the same conclusion in the case of adiabatic matter creation. This fact is in agreement with a 

previous result where we studied a full causal bulk viscous model through the renormalization 

group method arriving to the conclusion that in the long-time asymptotics the cosmological 

model tends to a flat ideal (non-viscous) Friedmann type geometry. 

We also notice that the x— parameter, i.e. the causal bulk viscous effect, weakly perturbs 

the FRW perfect fluid solution in such a way that the comoving entropy varies with t, 

while in the perfect fluid case, i.e. x = 0, the comoving entropy is constant. The same 

i 

considerations could be taken into account for the scale factor f(t) = f t :s <- 1 -^( u '+ 1 + >e ') , which 
is perturbed weakly for the bulk viscous parameter as we can see in the special case of 
7 = 1/2. The solutions obtained for the scale factor, energy density, entropy, etc. are similar 
"but different" to those obtained for a perfect fluid model where the x— parameter vanishes. 
We would like to point out that our model is thermodynamically consistent for the usual 
matter equations of state and valid for all 7— parameter, i.e. V7 G (0, 1) . In the same way, 
we can see that the viscous parameter helps us to get rid of the so-called entropy problem 
since in this model entropy varies with t. 

In the model considered, the evolution of the Universe starts from a singular state, with 
the energy density, bulk viscosity coefficient and cosmological constant tending to infinity. 
At the initial moment, t — 0, the relaxation time is r = 0. From the singular initial state, 
the Universe starts to expand, with the scale factor / as a 7-dependent function. The 
"constants" G, c and A are all decreasing functions of time. In the radiation predominance 
era we can take into account the effects of the bulk viscosity in all the quantities but when 
the universe enters a matter predominance era, these effects vanish and in this case our 
model behaves as a perfect fluid model with variable constants. The present model is not 
defined for s = 1, showing that in the important limit of small densities a different approach 
is necessary. 

The unicity of the solution is also proved by the investigation of the Lie group symmetries 
of the basic equation describing the time variation of the mass density of the Universe. 
Bulk viscosity is expected to play an important role in the early evolution of the Universe, 
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when also the dynamics of the gravitational and cosmological constants and the speed of 
light could be different. Hence the present model, despite its simplicity, can lead to a better 
understanding of the dynamics of our Universe in its first moments of existence. 
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